Nonlocal noise cross correlation mediated by entangled Majorana fermions 
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Due to their nonlocality, qubits nested in Majorana bound states may be the key to realize 
decoherence-free quantum computation. Majorana bound states could be achieved at the ends of 
a one-dimensional topological superconductor. However, when the bound states couple directly 
to electron reservoirs their nonlocal correlation is quenched by local Andreev reflections. Here we 
propose a scheme to generate nonlocal noise cross correlation between two well-separated quantum 
dots, mediated by a pair of Majorana bound states. Both positive and negative cross correlations 
can be obtained by tuning the gate voltages applied to the dots. Within a limited range of finite 
temperatures, the cross correlation is not suppressed by thermal fluctuations. Furthermore, we show 
how the local Andreev reflections suppress the noise cross correlation when multiple dot energy levels 
are coupled to the Majorana bound states. The measurable cross correlation is expected to serve as 
a sensitive indicator for the generation of Majorana fermions. 

PACS numbers: 03.75.Lm, 72.10.-d, 74.78. Na, 73.21.La 



I. INTRODUCTION 

Majorana fermions are exotic because they are their 
own antiparticlesii The search for Majorana fermions and 
controlling of Majorana bound states in solid-state de- 
vices are currently attracting increasing attention due to 
their potential application in quantum information^— 
Two well-separated Majorana bound states can define a 
nonlocal fcrmion level and its occupation encodes a qubit. 
This nonlocal topological qubit is robust against decoher- 
ence from local perturbations or interactions i^r— So far, 
a number of schemes have been proposed to realize and 
to detect the Majorana bound states*^— For instance, 
it was discovered that the Majorana bound states can be 
realized at the ends of a semiconductor nanowire with 
strong spin-orbit interaction in the proximity of an s- 
wave superconductor j 1 ^— Very recently, the signature 
for the formation of a spatially separated pair of Ma- 
jorana bound states have been observed experimentally 
in such a system, in terms of a zero-bias conductance 
peak i 21 ' 22 However, the existence of Majorana fermions 
is still controversial because the zero-bias conductance 
enhancement could also be the signature of diverse phe- 
nomena in mesoscopic physics . 23 ' 24 

A unique feature of Majorana bound states is that their 
superposition can form a nonlocal fermion level. Based 
on this nonlocality, the nonlocal current-current corre- 
lation in tunneling measurements would provide sup- 
porting evidence. Actually, the current noise cross cor- 
relation can reveal information related to particle fluc- 
tuations and is a powerful tool for studying different 
types of interactions and quantum statistics^ For free 
particles, fermions tend to induce negative noise cross 
correlation while bosons tend to induce positive one 
due to quantum statistics of indistinguishable identical 
particles^ In mesoscopic systems, positive cross correla- 
tion can be induced in the presence of Bardeen-Cooper- 
Schrieffer (BCS)-likc interactions^ - — dynamical chan- 
nel blockades ) 31 ' 32 and inelastic scattering— When a pair 



of Majorana bound states are coupled to two electron 
reservoirs directly, positive cross correlation may be in- 
duced by the nonlocal Andreev reflection, which means 
the injection of an electron into one bound state followed 
by the emission of a hole by the other— However, if 
two Majorana fermions are well separated, the nonlocal 
cross correlation might be absent ; 13 : 14 : 34 ' 35 which implies 
that it is hard to obtain the nonlocal signal of Majo- 
rana fermions by electrical measurements. This would 
become a serious obstacle for the application of Ma- 
jorana fermions in universal decoherence-free quantum 
computing. Furthermore, unlike Cooper pairs, in which 
the correlation is limited by the superconducting coher- 
ence length^2r— the correlation between two Majorana 
fermions is not restricted by their separation in space. In 
such a real nonlocal system, nonlocal noise cross correla- 
tion is expected to be much more important than local 
current autocorrelations such as shot noised Therefore, 
it is of fundamental interest to examine the noise cross 
correlation of two Majorana fermions in a long distance 
limit. 

In this paper, we propose a way to generate nonlo- 
cal current noise cross correlation mediated by a pair 
of Majorana bound states. Differing from a system 
where the Majorana bound states couple to electron 
reservoirs directl y 13 ' 14 ' 34 two quantum dots with single 
energy level are inserted between the Majorana bound 
states and reservoirs. In this case, local Andreev reflec- 
tions involving a single reservoir can be suppressed by 
weak dot-reservoir coupling, while crossed Andreev re- 
flections which split a Cooper pair over two reservoirs 
become dominant by strong coupling between the Majo- 
rana bound states and quantum dots. By tuning the gate 
voltages applied to the two dots, either positive or nega- 
tive nonlocal cross correlations can be induced in this hy- 
brid system of quantum dots and Majorana bound states. 
An experimental observation of the nonlocal noise corre- 
lation can serve as an alternative proof for the nonlocality 
of Majorana bound states. 
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The paper is organized as follows. In Sec. II, we pro- 
pose the model system for the double quantum dots and 
Majorana bound states and introduce the rate equation 
formulism employed to study the currents and their cor- 
relations. In Sec. Ill, we present the results for the noise 
cross correlation between the currents flowing through 
the double dots and discuss the underlying physical pro- 
cesses. We also consider the effects of thermal fluctuation 
and local Andrcev reflections enhanced by multiple dot 
levels. Finally, a summary is given in Sec. IV. 



II. MODEL AND FORMULISM 
A. Setup 

Our proposed setup is sketched in Fig. 1. A two- 
dimensional electron gas at a semiconductor hetero junc- 
tion, which is intensively studied in the context of both 
spin-orbit coupling and quantum dot, serves as the host 
for both the Majorana bound states and the dots. In 
the central region, spin-orbit coupling and strong Zee- 
man splitting give rise to a spinless energy band, which 
resembles ap + ip topological superconductor when adja- 
cent to an s-wavc superconductor wire due to the prox- 
imity effect ii2r— The nontrivial topological nature of the 
p + ip superconductor wire demands a pair of Majorana 
bound states to emerge at the two ends. Two quantum 
dots can be fabricated by depleting the electrons beneath 
the gate electrodes. Each quantum dot is coupled to a 
Majorana bound state and connected with the electron 
reservoir on its side. Only one spinless energy level is 
considered in each quantum dot by assuming sufficiently 
large level spacing and strong Zecman field. 




semiconductor heterojunction 



^ , Majorana . e/ 

dot1 J_ dot 2 \ kjf 

FIG. 1: Setup with the double quantum dots and Majorana 
bound states. An s-wave superconductor wire is deposited on 
the surface of a semiconductor heterojunction, and a magnetic 
field is applied along the axis of the wire. For a proper mag- 
netic field, two Majorana bound states appear at the ends 
of the topological superconductor. Two quantum dots are 
placed between the bound states and the reservoirs to block 
the local Andreev reflection between the Majorana bound 
states and electron reservoirs. The energy level of each dot 
can be modulated by applying a gate voltage. 



B. Hamiltonian 

The Hamiltonian for the system in Fig. [T] is given by 

t t * 

H = eidjdi + e 2 d 2 d 2 + -£a/?/i??2 

+(X* 1 d\ - Xnh)rn + i{X* 2 d\ + X 2 d 2 )r] 2 , (1) 

where Ej is the energy level in quantum dot i (i = 1,2) 
and di(d\) is the annihilation (creation) operator of elec- 
tron. The quantum dots arc coupled to the Majorana 
bound states with strength of Ai and X 2 . £a/ is the " cou- 
pling strength" between the Majorana fermions ??i(= r)\) 
and 172 (= 77!) in the Majorana bound states. The electron 
reservoirs and their coupling to the dots are described by 
the Hamiltonian 

h t = € hk4k c *k + J2(u4 k di + /i.e.), (2) 

1 A, 1 j A, 

where ti^k is the electron energy in the reservoir i and tj 
is the tunneling amplitude. 

It is helpful to replace the Majorana fermion operators 
by a fermion operator, 771 = f + f', T)i = i(f' — /), where 
/' creates a nonlocal fermion and P f = 0,1 counts 
the occupation of the corresponding statei 13 ' 17 In the 
new representation, the Hamiltonian in the central re- 
gion becomes^ 

H Q = € 1 d\d 1 +e 2 d\d 2 +e M {pf-^) 

+[\i(pd 1 + fdi) + X 2 (fU 2 - fd 2 ) + h.c.).(3) 

The Hamiltonian Ho can be solved in the space 
spanned by eight basis states \n\n 2 p), where n, = 0,1 
is the electron occupation number in the quantum dot i. 
For the Majorana bound states, the particle number in a 
superconducting states is not conserved, but the parity p 
can serve as the quantum number, with the values of even 
(e) or odd (o). For the double dots and Majorana bound 
states, the energy eigenvalues and eigenstates are given 
in two closed subspaces of definite parity. It is found that 
the states |00e), |10o), |01o), and | lie) form a closed even 
block, and other four states |00o), |10e), |01e), and llo) 
form a closed odd block. The eigenstates of even parity 
are 

\ei) = a1\00e) + 6f|10o) + cf\01o) + df |lle) 

and those of odd parity are 

\oi) = a°|00o) + 6°| 10e ) + c °|01e) + d°| n °) 

with I = 1,2,3,4 and where a,b,c,d are normalized 
superposition coefficients. For the two parities, their 
eigencquations are 



M o(e)^o(e) = E o(e)yo(e)^ 



(4) 
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In the absence of the electron tunneling between the 
quantum dots and reservoirs, there is no mixing between 
the states of different parities. The parity in the hy- 
brid system of Majorana bound states and dots can be 
varied by one-particle sequential tunneling between the 
quantum dots and the reservoirs. Thus, the even- and 
odd-parity sections are, thus, connected as 



(e l \d 1 \o v )=a?b° l ,+cfd° v , 



dl*c!„ 



(ei\d\\oi>) 

for the tunneling events through dot 1 and 

{e l \d 2 \o v )=ar4+brd° v , 
(e l \d\\o v ) = cra° l , + drb° l „ 



(7) 



(8) 



for dot 2. The dot-reservoir tunneling does not mix 
the states of the same parity, i.e., (e;|dL 2 Jej') = and 

(oz|4 (2) |of) = 0. 



C. Rate equations 

We study the tunneling currents and their noise cross 
correlation in the sequential tunneling regime with the 
help of the rate equation method^ Formally, there are 
two types of approaches: the generalized master equa- 
tions and diagonalized rate equations. In the generalized 
master equations, the density matrix contains not only 
the population of each state (may not be the eigenstates) 
but also the coherent terms between any two states. The 
diagonalized rate equation approach is to, first, diago- 
nalize the Hamiltonian of central region and then write 
down the rate equations for the population of central 
eigenstates. For weak coupling between the dots and 
reservoirs, it is safe to construct the density matrix p 
with the eigenstates of the hybrid system of Majorana 
bound states and dots then write down the rate equa- 
tions for the population of the system eigenstates. It has 
been shown that when the inner interaction of the cen- 
tral region is much stronger than the dot-lead coupling 
strength, the results obtained from two approaches are in 



close agreement with each otheri 37 ! 38 Therefore, we take 
A.; 3> U in the calculation to ensure that the diagonal- 
ized rate equations work well. The time evolution of the 
density matrix p(t) is given by the rate equations 



dp(t)/dt = Wp{t), 
where the rate matrix elements ar o 37 ' 39 

W in = J2 T ilf( A i'i + ^)\(l'\di\l)\ 2 

i 

+f(A n -f H )\(l'\dl\l}\ 2 ] 

for I ^ I' and 

N 



(9) 



(10) 



(11) 



with l,V g {|oj),|ei}}. Here f(u) = [1 + e UJ ' kBT ]- 1 is 
the Fermi distribution function, A;/; = Ey — Ei is the 
Bohr frequency of the transition from |Z) to |/'), Ei is 
the eigenencrgy of \l), and pi is the chemical potential 
in reservoir i. In the wide-band limit approximation, the 
coupling between the i-th dot level and its reservoir is 
denoted by Ti = 2ir\ti\ 2 Vi with Vi the spinless density of 
states near the Fermi surface of reservoir i. The current 
Ii flowing through dot i is calculated by 



(12) 



where T" 1 is the matrix form of the current operator and 
its elements are given by 



n 



VI 



r, [f(A in 



f iM)\{l'\di 



(13) 



The noise cross correlation of two tunneling currents 
can be studied by the Fourier transform of the current- 
current correlation function 20 



S 12 {oj) 



dte- t ((7 1 (t)7 2 (0))-(/ 1 )(7 2 )) 1 (14) 



where (• • • ) represents the thermodynamic average. Fur- 
thermore, the current-current correlation function of the 
currents I\ and 7 2 in the u> space can be expressed as 



(h{t)I 2 {0)) u = [f if ( W )f 2 p (0) + f 2 f (-W)f XP 



(0) 



(15) 



where T(±w) = (=pio;I — W) and I is the unit matrix. 

In the following discussion, we focus on the situation 
that the two Majorana fermions are well separated so 
that cm = 0. The two dots are assumed symmetrically 
coupled to the reservoirs Ti = Tq. Considering the sym- 
metric couplings between Majorana bound states and 
dots on the opposite ends, we take A; = Aq. 
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III. RESULTS AND DISCUSSIONS 

A. Tunable nonlocal noise cross correlation 

The zero-frequency noise cross correlation £12 is pre- 
sented in Fig. 2 as a function of dot levels t\ and £2, 
where the voltage V = —fa is symmetrically applied to 
both dots. The correlation shows a quadrupole pattern of 
four regions, two positive and two negative. At the point 
ei = £2 = 0, the cross correlation is zero. Away from this 
point, the cross correlation with tunable signs can be in- 
duced by varying the dot energy levels. Along the diago- 
nal directions, the correlation reaches the maxima. When 
both quantum dot levels are higher or lower than the 
chemical potential (set as zero energy point) in the p + ip 
superconductor, the correlation is negative. Otherwise, 



the correlation is positive. Similar patterns were also ob- 
tained in the capacitively coupled double dots^ where 
the cross correlation arises from the competition between 
the inter- and intradot Coulomb interactions. Here the 
sign reversal of the cross correlation is mediated by the 
Majorana bound states. Competition and coordination 
between different tunneling channels are modulated by 
the occupation of Majorana bound states. As shown in 
Fig. 2(a), two negative cross-correlation branches lie in 
the competition region asymmetrically and the maximum 
negative S12 appears at the point = fa. Away from 
this point, electron tunneling into (out of) the dots be- 
comes easier when ei < [>)fa. This weakens the com- 
petition between tunneling currents on the opposite ends 
and results in the decrease of noise cross correlation. 
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FIG. 2: Tunable noise cross correlation, (a) Zero-frequency current noise cross correlation as a function of two dot energy 
levels ei and £2- The parameters are taken as em = 0.0, ksT = 2.0, To — 0.2, Ao = 5.0, and fa — —2.0. (b) When £1,62 < 0, 
the system tends to reside in the state \n1n2p) = |Hp), where m £ {0, 1} is the occupation in dot i = 1,2, the parity p of the 
Majorana bound states could be even or odd. | lip) could evolve into either |Qlp) or |10p), where p stands for the opposite parity 
to p. For both |01p) and |10p), electron tends to tunnel in one dot and out of the other. These opposite tunneling processes 
lead to the negative cross correlation. Similarly, the correlation is also negative in (e) when £i,£2 > 0. (c) When e\ < < £2, 
the possible transitions are |10p) «-»• |00p), |llp), which tend to synchronize the tunnelings in two dot-reservoir channels and 
give rise to the positive correlation. Similar processes also give the positive cross correlation in (d) when £2 < < £1. 



B. Physical processes 

To analyze the physical processes beneath the correla- 
tion, we return to the /-representation in Eq. (3) , which 
shows three types of transitions between the Majorana 
bound states and quantum dots: pairing creation pd\, 
pairing annihilation fdi, and direct tunneling f'di and 
d\f. We, first, take Fig. 2(b), for example. By tun- 
ing the gate voltages, one can lower the energy levels 



in both dots than the chemical potentials in the reser- 
voirs (ei <C fa). Consequently, the two dots are occupied 
most of the time and the system is in the state | lip) 
(where p can be even e or odd o). In this state, ther- 
mal fluctuations can induce tunneling currents in both 
dot-reservoir channels, but the currents are uncorrelatcd 
when they are not related to the Majorana bound states. 
For ei -c fa, the state population of \00p) is rather small 
and almost does not participate in the state evolution. 
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For weak dot-reservoir coupling, the states can evolutc 
many times before parity is changed by an electron tun- 
neling into or out of the reservoirs. The state |llp) can 
evolve into |10p) or \01p) by direct tunneling if p = e 
or pairing annihilation if p = o, where p = e if p = o 
and vice versa. In either case, the occupations in the two 
dots change from |11) into 1 10) or |01). As a result of 
the interaction between the dots and Majorana bound 
states, electrons tend to tunnel into one dot or out of the 
other. These opposite tunneling behaviors lead to the 
negative cross correlation between the currents of two 
dot-reservoir channels. Similarly, for the case q 3> /Ltj 
[Fig. 2(e)], the initial state |00p) can evolve into |10p) 
or |01p), which also gives rise to the negative cross cor- 
relation. In contrast, when the zero energy point in the 
region of Majorana bound states is much higher than 
one dot energy level while much lower than the other, 
the possible transitions become from | lOp) to |00p) or 
| lip), as shown in Fig. 2(c), and from Olp) to |00p) or 
| lip), as shown in Fig. 2(d). These transitions tend to 
synchronize the electron tunneling processes in both dots 
and make the cross correlation positive. 

It is indicated in Fig. 2(a) that the noise cross correla- 
tion S12 reaches the maximum values but with opposite 
signs on the two diagonal lines e\ = ti and ei = —62. In 
these two cases, we can understand the nonlocal correla- 
tion from analytical results. In the case e\ = £2 = cq, the 
eigenvalues of Hq can be found as 

El± = \ U ± V^o + SA^) , (16) 
and the corresponding eigenstates are 



*? ± = .|io P -) + |oip-) + £o± V° + 8Ao |iip), 



*S± = £oT V" ° |oo P -)-|io P -) + .|oi P -), (17) 

ZAq 

where q — 1 for even parity and <r = — 1 for odd parity. 
When a relatively small voltage is applied to both dots, 
the eigenstate (^2-^) ^ s ^ ne ground state in the re- 

gion of Majorana bound states and dots if the dot energy 
levels are lower (higher) than the chemical potential in 
the nanowire, i.e., eo < (>)0. For eo < 0, the stationary 
population of ^>°^ is dominant and the transition be- 
tween these two degenerate states could be caused by the 
tunneling between dots and reservoirs. When there is no 
interaction between the Majorana bound states and dots, 
the ground state is the one that the each dot is always 
singly occupied most of the time. The dot-Majorana in- 
teraction correlates the two dots nonlocally and induces 
the transition events |llp) -H- Olp) or lOp) . This means 
that the role of the coupling between the Majorana bound 
states and dots is to turn the same state in the two dots 
into opposite simultaneously. Therefore, the tunneling 



events between two channels tend to repulse each other, 
which gives rise to a negative noise cross correlation. Sim- 
ilarly, the ground state of the system is $2- m tne case 
of eo > 0. The interaction between the Majorana bound 
states and dots causes the transition events |00p) f-> |01p) 
or 1 10p) . The noise cross correlation is also negative in 
this case. 

When the two dot levels have opposite energies ei = 
— e 2 = £o, the eigenvalues of Hq are 

^ ± = ^o±V^ + 8A^, (18) 
and the corresponding eigenstates are 



«S± = IQOP) + £ ° V° + aA ° |10p)+^|llp). (19) 

AAq 

It is clearly shown that the coupling between the Ma- 
jorana bound states and dots turns the opposite states 
of double dots into the same simultaneously, e.g., Olp) 
(|10p)) f-» |00p) or I lip). The presence of Majorana-dot 
coupling tends to correlate the tunneling events in both 
dots in a bunching way. Therefore, the noise cross corre- 
lation is always positive in this case. 

C. Effects of thermal fluctuation 

In the sequential tunneling regime, the effect of ther- 
mal fluctuation is reflected in the Fermi distribution of 
tunneling electrons in both leads. The reservoirs are not 
affected by the coupling to the dots so they remain in 
their respective thermal equilibrium. With the decrease 
of temperature, fewer electrons away from the Fermi en- 
ergy contribute to the current and thermal fluctuation- 
induced tunneling is suppressed. Figure 3 shows that the 
noise cross correlation S12 exhibits quite different pat- 
terns for different temperatures ksT. At low tempera- 
tures, the decay time of the lead correlation is mainly 
determined by the bias voltage, while the contribution 
from the quasiparticle lifetime is irrelevant] 40 ' 41 In the 
calculation, we take a bias voltage much larger than the 
dot-reservoir strength, which ensures that the results at 
relatively low temperatures are still justified. Figure 3(a) 
represents the results when the temperature is lower than 
the dot-lead coupling strength. In the region of negative 
values, the cross correlation decreases to zero for e% > 
and shrinks to an elliptical region centering at eo = 
Considered individually, the Majorana-dot Hamiltonian 
Hq has the particle-hole symmetry. However, when the 
reservoirs with tunable chemical potentials are connected 
to the dot, the particle- hole symmetry is broken. In 
the calculation, we take the chemical potential in both 
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FIG. 3: The current noise cross correlation as a function 
of two dot energy levels ei and £2 for different temperatures 
k B T. (a) k B T = 0.5, (b) k B T = 2.0, (c) k B T = 5.0, and (d) 
k B T = 10.0, The parameters are taken as To = 1.0, Ao = 5.0, 
and Hi = —2.0. 



leads [i = —2.0 to generate two currents. Here the bro- 
ken particle-hole symmetry is attributed to the finite 
bias voltage and the cross correlations thus differ in the 
regimes of e,; > and < 0. With the increase of tem- 
perature, four regions in the cross-correlation pattern be- 
come wider and more symmetric about the point (0,0). 
From Eq. (fTU)) it is found that the tunneling probability 
of electrons is governed by the Fermi distribution func- 
tion /(A// j ±/Uj). When fc^T is much larger than Ao, the 
tunneling probability is determined by the ratio ei/ksT. 
Therefore, the area of cross correlation is proportional 
to the energy scale of the temperature approximately, as 
shown in Fig. 3. 

It is interesting that the temperature does not suppress 
the strength of cross correlation, which differs from the 
case of capacitively coupled double dots^ For interdot 
Coulomb interaction-induced cross correlations in dou- 
ble dots, each dot is connected to two electron reservoirs. 
At a relatively high temperature, the thermal fluctuation 
helps the electrons tunnel into or out of the dots more 
easily. The (anti-)blockade effect between two tunneling 
channels is, thus, weakened, leading to a suppression of 
noise cross correlation. However, for the present device, 
one of the reservoirs is replaced by one end of the topolog- 
ical superconductor. The zero-mode Majorana fcrmions 
serve as the source or drain, in which only one fermionic 
level takes part in the transport. Although the temper- 
ature determines the state population of the system, the 
mechanism of Majorana mediated cross correlation is in- 
dependent of the temperature. 



D. Suppression of nonlocal noise correlation by 
enhanced local Andreev reflection via multiple dot 
levels 

In the setup where the Majorana bound states cou- 
ple directly to the reservoirs, the noise cross correlation 
disappears (S12 = 0) for vanishing eM ' 13 ' 14 ' 34 Strictly 
speaking, cm must be according to the definition that 
the Majorana fermions are their own antiparticles. In the 
present device the direct couplings between the Majorana 
bound states and reservoirs are blocked by the quantum 
dots. The weak dot-reservoir coupling suppresses the lo- 
cal Andreev reflection between the reservoirs and Majo- 
rana bound states and enhances the nonlocal cross cor- 
relation. To test this picture, we also study two cases in 
which multiple dot energy levels are coupled to the Majo- 
rana bound states: (1) only one end of the topological su- 
perconductor is coupled to multiple levels and; (2) both 
ends are coupled to multiple levels. The two cases are 
schematically shown in Figs. 4(a) and 4(b), respectively. 
The energy levels provide more conducting channels be- 
tween the reservoirs and Majorana bound states, thus 
enhancing the local Andreev reflection. As the number 
N of the energy levels approaches infinity, the model be- 
comes equivalent to the case when the Majorana bound 
states are coupled to electron reservoirs directly. As a 
result, the current noise cross correlation is expected to 
vanish. 



(a) 



N 



(b) 



k A 



N 



FIG. 4: Schematic view for the cases when multiple levels 
are coupled to the Majorana bound states at the ends of the 
topological superconductor. Multiple levels are coupled to (a) 
only one end of the topological superconductor and (b) both 
ends. Each level is assumed to have the same energy and 
coupling to the Majorana bound states. 



Figure 5 demonstrates the noise cross correlation 
as a function of the energy level ed for different level num- 
bers iV. For simplicity, it is assumed that the levels are 
of the same energy and coupling strength to the Majo- 
rana bound states. For comparison, Fig. 5(a) presents 
the result for the first case in Fig. 4(a), in which only 
one end couples to multiple levels, while Fig. 5(b) gives 
the result of the second case in Fig. 4(b) where both 
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FIG. 5: Current noise cross correlation as a function of energy 
level td in the presence of multiple levels, where N counts 
the level number. The N levels couple to only one end (a) 
and both ends (b) of the topological superconductor, where 
the coupling strength between the N levels and reservoirs is 
set as r = To/TV. The parameters are taken as em = 0.0, 
k B T = 2.0, T = 0.2, A = 10.0, and fa = -2.0. 



ends couple to multiple levels. For N = 1, the model 
reduces to the case along t\ = e 2 in Fig. 2. In this case, 
the maximum negative S12 can be achieved. It is found 
that with increasing number of energy levels, the cross 
correlation S12 decreases rapidly in both cases. In the 



presence of multiple levels, the electron tunneling caused 
by Andreev reflection does not have to always involve 
both electron reservoirs but can occur via two channels 
in the same reservoir. In this way, the cross correlation 
between the currents flowing through two reservoirs is 
suppressed. Since the Andreev reflection can occur via 
any two levels, it is more likely to occur at the end of 
more energy levels for the case in Fig. 4(a). As a result, 
the suppression of S12 by extra levels in Fig. 5(a) is much 
stronger than that in Fig. 5(b) as the level number N 
increases. 

IV. SUMMARY 

We propose a setup to generate nonlocal noise cross 
correlation between the currents flowing via two well- 
separated quantum dots, mediated by a pair of entan- 
gled Majorana fermions. Each dot is inserted between a 
Majorana bound state and a nearby electron reservoir to 
suppress the local Andreev reflection between them. In 
this case, the crossed Andreev reflection over two dots be- 
comes dominant due to the strong coupling between the 
Majorana bound states and dots. It is demonstrated that 
the nonlocal cross correlation with tunable signs can be 
induced by the Majorana-dot coupling and modulated by 
the dot energy levels. Because the nonlocal correlation 
cannot be induced by other mechanisms, these measur- 
able noise cross correlations can serve as a unique signa- 
ture for the formation of Majorana bound states. It is 
found that the cross correlation is not suppressed with 
the increase of temperature and the area of correlation 
is proportional to the temperature approximately. More- 
over, we study the case when multiple dot energy levels 
are coupled to the Majorana bound states. It is shown 
that with the increase of energy level number, the local 
Andreev reflection becomes enhanced and suppresses the 
crossed Andreev reflection considerably. As a result, the 
current noise cross correlation vanishes gradually. 

We thank T. K. Ng and K. T. Law for insightful dis- 
cussions. This project was supported by the Research 
Grant Council under Grant No. HKUST3/CRF/09, the 
Foundation for Innovative Research Groups of the NSFC 
under Grant No. 61021061, and the NSFC (Grants No. 
11004022 and No. 61006081). 



1 E. Majorana, Nuovo Cimento 14, 171 (1937). 

2 A. Stern, Nature (London) 464, 187 (2010). 

3 C. Nayak, S. H. Simon, A. Stern, M. Freedman, S. Das 
Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

4 G. Moore and N. Read Nucl. Phys. B 360, 362 (1991). 

5 A. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 

6 L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 

7 J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 
Phys. Rev. Lett. 104, 040502 (2010). 

8 J. Alicea, Y. Oreg, G. Refael, F. von Oppen, and M. P. A. 
Fisher, Nat. Phys. 7, 412 (2011). 

9 J. Linder, Y. Tanaka, T. Yokoyama, A. Sudb0, and N. 



Nagaosa, Phys. Rev. Lett. 104, 067001 (2010). 

10 Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 
105, 177002 (2010). 

11 R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. 
Lett. 105, 077001 (2010). 

12 J. Alicea, Phys. Rev. B 81, 125318 (2010). 

13 C. J. Bolech and E. Dernier, Phys. Rev. Lett. 98, 237002 
(2007). 

14 K. T. Law, P. A. Lee, and T. K. Ng, Phys. Rev. Lett. 103, 
237001 (2009). 

15 D. E. Liu and H. U. Baranger, Phys. Rev. B 84, 201308(R) 
(2011). 



8 



16 A. Golub, I. Kuzmenko, and Y. Avishai, Phys. Rev. Lett. 
107, 176802 (2011). 

17 M. Leijnse and K. Flensberg, Phys. Rev. B 84, 140501 (R) 

(2011) . 

18 M. Wimmer, A. R. Akhmerov, J. P. Dahlhaus, and C. W. 
J. Beenakker, New J. Phys. 13, 053016 (2011). 

19 K. Flensberg, Phys. Rev. Lett. 106, 090503 (2011). 

20 L. Fu, Phys. Rev. Lett. 104, 056402 (2010). 

21 V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. 
M. Bakkers, and L. P. Kouwenhoven, Science 336, 1003 

(2012) . 

22 M. T. Deng, C. L. Yu, G. Y. Hua ng, M. Larsson, P. Caroff, 
and H. Q. Xu. larXiv:1204.4130l 

23 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. 
Abusch-Magder, U. Meirav, and M. A. Kastner, Nature 
(London) 391, 156 (1998). 

24 S. M. Cronenwett, H. J. Lynch, D. Goldhaber-Gordon, 
L. P. Kouwenhoven, C. M. Marcus, K. Hirose, N. S. 
Wingreen, and V. Umansky, Phys. Rev. Lett. 88, 226805 
(2002). 

25 Y. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 

26 W. D. Oliver, J. Kim, R. C. Liu, and Y. Yamamoto, Sci- 
ence 284, 299 (1999). 

27 L. Hofstetter, S. Csonka, J. Nygard, and C. Schonenberger, 
Nature (London) 461, 960 (2009). 

28 J. Wei and V. Chandrasekhar, Nat. Phys. 6, 494 (2010). 

29 A. Freyn, B. Dougot, D. Feinberg, and R. Melin, Phys. 



Rev. Lett. 106, 257005 (2011). 

30 P. Recher, E. V. Sukhorukov, and D. Loss, Phys. Rev. B 
63, 165314 (2001). 

31 A. Cottet, W. Belzig, and C. Bruder, Phys. Rev. Lett. 92, 
206801 (2004). 

32 D. T. McClure, L. DiCarlo, Y. Zhang, H. A. Engel, C. 
M. Marcus, M. P. Hanson, and A. C. Gossard, Phys. Rev. 
Lett. 98, 056801 (2007). 

33 C. Texier and M. Biittiker, Phys. Rev. B 62, 7454 (2000). 

34 J. Nilsson, A. R. Akhmerov, and C. W. J. Beenakker, Phys. 
Rev. Lett. 101, 120403 (2008). 

35 G. Striibi, W. Belzig, M. S. Choi, and C. Bruder, Phys. 
Rev. Lett. 107, 136403 (2011). 

36 S. Tewari, C. Zhang, S. Das Sarma, C. Nayak, and D. H. 
Lee, Phys. Rev. Lett. 100, 027001 (2008). 

37 C. Poltl, C. Emary, and T. Brandes, Phys. Rev. B 80, 
115313 (2009); H. Z. Lu, B. Zhou, and S. Q. Shen, Phys. 
Rev. B 79, 174419 (2009). 

38 C. Emary, C. Poltl, and T. Brandes, Phys. Rev. B 80, 
235321 (2009). 

39 H. Bruus and K. Flensberg, Many-Body Quantum Theory 
in Condensed Matter Physics (Oxford University Press, 
Oxford, UK, 2004). 

40 A. Mitra and A. J. Millis, Phys. Rev. B 76, 085342 (2007). 

41 C. Timm, Phys. Rev. B 77, 195416 (2008). 



